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ABSTRACT 

The afterglow of a cosmological Gamma-Ray Burst (GRB) should appear on the 
sky as a narrow emission ring of radius ~ 3 x 10^^ cm (t/day)^/^ which expands faster 
than light. After a day, the ring radius is comparable to the Einstein radius of a 
solar mass lens at a cosmological distance. Thus, microlensing by an intervening star 
can modify significantly the lightcurve and polarization signal from a GRB afterglow. 
We show that the achromatic amplification signal of the afterglow flux can be used 
to determine the impact parameter and expansion rate of the source in units of the 
Einstein radius of the lens, and probe the superluminal nature of the expansion. If the 
synchrotron emission from the afterglow photosphere originates from a set of coherent 
magnetic field patches, microlensing would induce polarization variability due to the 
transient magnification of the patches behind the lens. The microlensing interpretation 
of the flux and polarization data can be confirmed by a parallax experiment which 
would probe the amplification peak at different times. The fraction of microlensed 
afterglows can be used to calibrate the density parameter of stellar-mass objects in the 
Universe. 

Subject headings: cosmology: theory - gamma rays: bursts - gravitational lensing 

1. Introduction 

The recent discovery of delayed X-ray (van Paradijs et al. 1997), optical (Bond 1997; 
Djorgovski et al. 1997; Mignoli et al. 1997), and radio (Frail et al. 1997) emission, over hours 
to several months following 7-ray bursts (GRB) established a new class of variable sources in 
astronomy. Of particular significance is the detection of Fell and Mgll absorption lines at a 
redshift of z = 0.835 in the optical spectrum of GRB970508 (Metzger et al. 1997), which confirmed 
the extragalactic origin of this burst. Since the source redshift must be higher than the absorber 
redshift, its required optical luminosity exceeds that of a supernova by several orders of magnitude. 
Thus, GRB afterglows might be detectable out to high redshifts. One could then use the signatures 
of absorption in the optical band (Metzger et al. 1997; Djorgovski et al. 1997), scintillations in 
the radio regime (Goodman 1997), or gravitational lensing (Gould 1992; Mao 1993) by intervening 
material along the line-of-sight, to study the intrinsic properties of afterglow sources. 
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Afterglows are most naturally explained by models in which the bursts axe produced by 

rclativistically expanding fireballs (Paczynski & Rhoads 1993; Meszaros &: Rees 1997; Vietri 1996; 
Waxman 1997a, b; Wijers, Rees, & Meszaros 1997; Vietri 1997; Sari 1997). On encountering an 
external medium, the relativistic shell which emitted the initial GRB decelerates and converts its 
bulk kinetic energy to synchrotron radiation, giving rise to the afterglow. The combined radio 
and optical data imply that the fireball energy is ~ io5i-52 gj.g j^^g relativistic beaming, 
the emission region seen by an external observer, occupies an angle ~ I/7 relative to the center 
of the explosion, where 7 is the Lorentz factor. This region appears to expand faster than the 
speed of light and occupies an angle of ~ 0.1 — 10^ micro-arcscconds on the sky (or a physical 
size of ~ 10^^-10^^ cm). Due to the smallness of this angular size, it is difficult to resolve the 
afterglow source by terrestrial telescopes. However, the lensing zone of a solar mass lens located 
at cosmological distances occupies a micro-arcsecond on the sky (hence the term "microlensing" ) , 
and thus offers the unique opportunity for resolving GRB sources during their afterglow phase. 
Because of the superluminal expansion of the source, any (non-relativistic) peculiar velocity of the 
lens relative to the source can be ignored. The amplification peak of a microlensing event lasts 
for only ^ day, after which the net amplification weakens as the source size grows larger than the 
Einstein radius of the lens. The short duration of a microlensing event could therefore provide 
a test for the high Lorentz factor of the afterglow photosphere, which is predicted by all fireball 
models (for comparison, the variations due to peculiar velocities in microlensing events of steady 
sources take decades rather than days). 

The rapid expansion and deceleration of the fireball causes a sharp decline in its surface 
brightness as a function of time. Since emission along the line-of-sight to the source center suffers 
from the shortest geometric time-delay, it occurs at larger radii and appears dimmer relative to 
slightly off'-axis emission. At any given time, the source is expected to appear as a narrow ring 
of radius R/'~f and a width of order a tenth of this radius (Waxman 1997c). The outer cut-off 
is set by the sharp decline in relativistic beaming outside the ring. As the ring crosses a lens, 
its magnification adds a sharp peak to the otherwise smooth light curve of the afterglow. The 
sharpness of the peak depends on the thickness of the radiating gas layer behind the shock and on 
the shock deceleration rate. Microlensing could therefore provide important information about the 
structure and dynamics of the afterglow photosphere. 

The probability for stellar microlensing of a source at a redshift Zj, ~ 1 is ~ 0.1il^6^ (Press &: 
Gunn 1973; Gould 1995), where is the mean density of stellar-mass objects in the Universe in 
units of the critical density, and h is the impact parameter of the source relative to the lens in 
units of the Einstein radius. The known population of luminous stars amounts to Jl,^ ~ 5 x 10~^ 
(Woods & Loeb 1997), and implies that most cosmological sources are separated from stellar 
lenses by 6 40. The typical impact parameter is smaller by an order of magnitude if the dark 
matter is made of Massive Compact Halo Objects (MACHOs) as Galactic microlensing searches 
suggest (Alcock et al. 1996). 

In this paper we examine the question whether a stellar mass lens can resolve the predicted 
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properties of afterglow photospheres. For concreteness, we derive numerical results for the fireball 
emission model of Waxman (1997b, c). Since the afterglow occurs long after the explosive energy 
release, its properties are not sensitive to the spatial or temporal details of the point explosion that 
triggered the GRB. However, our adopted emission model is by no means a unique interpretation 
of the existing afterglow data (see, e.g. Vietri 1997 or Paczyhski 1997); in fact, a future detection 
of a microlensing signal could serve to discriminate among competing afterglow models. 

In §2 we describe our model for GRB afterglows and characterize both the intensity and 
polarization signals that would result from a microlensing event. The numerical results and their 
implications are discussed in §3. Finally, §4 summarizes the main conclusions from this work. 



2. Source Model and Microlensing Signatures 
2.1. Source Model 

To illustrate the effects of microlensing on GRB afterglows we need to specify the evolution of 
the source size and spectral intensity with time. We adopt the scaling laws for the expansion and 
emission of a relativistic fireball which decelerates in a uniform ambient medium (Waxman 1997b). 

In the fireball model, a compact (~ 10^~^ cm) source, releases an energy of E ^ 10^^ ergs 
over T ^ 10^ seconds with a negligible baryonic contamination 10"^ Mq). The high energy 
density at the source results in an optically thick pair plasma that expands and accelerates to 
relativistic velocities. After an initial acceleration phase, the thermal energy is converted to kinetic 
energy of the protons. A cold shell of thickness cT is formed and continues to expand. Internal 
shell collisions as a result of unsteady source activity could convert part of kinetic energy into 
radiation and yield the primary GRB emission via synchrotron emission and inverse-Compton 
scattering (Paczyriski & Xu 1993; Meszaros & Rees 1994; Sari & Piran 1997). As the cold shell 
expands, it impacts on the surrounding medium and drives a relativistic shock also in it; this shock 
continuously heats fresh gas and accelerates relativistic electrons which produce via synchrotron 
emission the delayed radiation observed on time scales of hours to months. Following Waxman 
(1997b), the radius of the shock at observed time t is given by 

R{t) « 8.7 X 10^^ Eli^-^^hll'^ cm (1) 

while its Lorentz factor is 

7W = /^-21eVW/«. (2) 

Here £^52 is the fireball energy in units of 10^^ ergs, ni is the ambient gas density in cm~^ and thr 
is the observed time in hours. As mentioned in §1, most of the emission is seen from a narrow ring 
of radius 

P^y-> ~ TJ+\ ^ ^52 "-1 *hr cm . (3j 
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The width of the ring is a fraction W ~ 10% of its radius ps if the thickness of the radiating layer 
behind the shock is determined by the shock hydrodynamics in a self-similar expansion (Waxman 
1997c). A thicker radiating layer (e.g. due to a large gyroradius of the radiating electrons) would 
result in a wider ring. In §3, we will show numerical results for different choices of W. These 
expressions are valid also for a jet geometry as long as the opening angle of the jet is ^ I/7. 

The X-ray, optical and radio emission following the 7-ray burst can be modelled as 
synchrotron emission from a power-law population of electrons within the heated shell behind the 
expanding shock. Under the assumption that the magnetic field energy density in the shell rest 
frame is a fraction S,b of the equipartition value, and that the power-law electrons carry a fraction 
of the dissipated energy, the observed frequency at which the synchrotron spectral intensity of 
the electrons peaks is 

where Zg is the cosmological redshift of the source. The observed intensity at i/^ is 

^i'{^Y'E,,mJy. (5) 

If the distribution of electron Lorentz factors follows a power-law, dNg/dje oc 7^7^, with a low 
energy cut-off set by ^e, then the observed intensity as a function of frequency, i^, obeys, 

F^it) = F,^ [v/v^it)]-'^ , (6) 

where Vm is the emission frequency of the electrons at the low-energy cut-off. The variation in u^a 
across the finite width of the ring can be ignored for W <^\. The typical parameter values which 
are required to fit the afterglow data are ~ 0.1, ~ 0.1, and p ~ 2, so that a ~ 0.5 for v > Vm 
and a = —1/3 for u < Vm- 
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2.2. Flux Amplification Due to Microlensing 

We now consider a point lens of mass M and redshift zi which happens to be located near the 
line-of-sight to an expanding fireball. We denote by r] the impact parameter of the source center 
relative to the observer- lens axis. For simplicity we assume that the source has a uniform surface 
brightness in a ring of radius Ps{t) [given by equation dH)] and a width W ps{t). The flux seen by 
the observer is 

F'r'%Rs{t),W,b] = F^{t) ,ji[Rs{t),WM , (7) 

where Rg = Ps/^E, and is the Einstein radius of the lens projected on the source plane, 

te = \/ {AGM/c^){DsD\^/ D\), with Di,Ds and Dig being the angular diameter distance to the lens, 

to the source, and from the lens to the source, respectively. These distances all depend on the 
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cosmological parameters. In this paper we assume Q = 1, A = 0, and Hq = 50 km s ^ Mpc ^. 
The magnification factor for a normahzed lens-source separation 6 = r^/rs is, 

where ^{Rs,b) is the magnification for a uniform disk of radius Rs (Schneider, Falco, & Ehlers 
1992), 



nRl 



/ dr^^arccos^— 1 + HiR^ - b)-{Rs - b)J {R^ - b)^ + i . 

J\b-Rs\ 2r6 2 \' 

(9) 

Here H{x) is the Heaviside step function. The integral in equation (^) can be expressed more 
explicitly as a sum of elliptic integrals (Witt & Mao 1994). Other analytic results exist for more 
general surface brightness distributions (Heyrovsky & Loeb 1997). 



2.3. Polarization Variability Due to Microlensing 

If the afterglow photosphere contains a finite set of discrete patches, each having a coherent 
magnetic field distribution, then the emergent synchrotron radiation will be polarized. For a 
power-law distribution of electron energies with an index p, the degree of linear polarization in 
each coherent patch is given by (Rybicki & Lightman 1979) 

p + 7/3 

For the inferred value of p ~ 2 (Waxman 1997b), 11 ~ 0.7. A microlens capable of resolving the 
source, could then provide useful information about its magnetic field structure. 

To illustrate the effect of microlensing on polarization we adopt a toy model in which the 
emission ring is divided into a set of independent segments, each having a coherent distribution of 
magnetic field lines. The polarization in each segment is then modelled as a traceless symmetric 
2x2 tensor with a random orientation angle, and a contraction {PapP^^Y^"^ = IIj given by 
equation (|lO|). To simplify the computation, we subdivide the emission ring into segments 
of equal area and nearly square shape. To each segment we assign a randomly oriented linear 
polarization. For the sake of concreteness, we assume that the number of segments and the 
orientation of their polarization stays constant during the lensing event. This assumption is 
reasonable since the effect of lensing peaks during the short period of time when the ring crosses 
the lens (which is smaller than the expansion time by a factor W <^1). 

The net observed polarization is then given by 



{P) = ^^^^ , (11) 
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where Ai and Pi are the area and polarization tensor of the i-th segment, and 

p ^ Il_ / cos2(t)i sin2(/>j \ 
' ^ V sin 2(pi - cos 2(j)iJ ' ^ ^ 

with a random orientation angle, < (pi < 27r. 

We first consider the case where there is no lensing, in which Ai = Aq = const for i = 1, ...A^. 
The two components of the net polarization are then given by 



EiliH sin 20,^0 n ^ 



{P),y = {P)y, = ^^=1 — -v^^--u = ^ysin2<p, . 



(13) 



Clearly, the resulting polarization (P) = y 2((P)2^ + {P)%y) approaches zero for large and is 
time independent. 

Let us now consider the situation where a lens is located at a projected position {xi,yi) with 
respect to the center of the source, so that b = {xf + yfy^'^. The observed polarization is still given 
by equation (11) but due to the stretching caused by lensing, the areas {Ai}^L-^ of the various 
segments are no longer equal, 

where (C, 0) are polar coordinates centered on the source, and the integrand is the point-source 
amplification factor at an impact parameter d = [(^cos^ — x;)^ + (Csin^ — yi)'^]^^'^. The integral 
is taken over the unlensed area of the segments. Because the size and position of the various 
segments relative to the lens change with time, the observed polarization will vary during a 
microlensing event. 



3. Numerical Results 

3.1. Flux Amplification 

The solid lines in Figure 1 show the unlensed lO^^Hz flux of an afterglow according to 
equation (^) and the parameter choices mentioned below that equation. The broken lines show 
the effect of microlensing on the observed flux, according to equation (|7|), for different choices of 
the ring's fractional width W and impact parameter b. 

The qualitative features of the microlensing signature on the afterglow lightcurve are as 
follows: 



-7- 



(a) All wavelengths show the same amplification profile as a function of timeQ. While the 
amplification peak occm's on the rising side of the lightcurve in the radio, it appears on its 
declining side in X-rays, and might show on both sides of the break in the optical [see Eq. (Q) 
for the timing of the peak at a given frequency]. The larger b is, the easier it becomes to 
detect the amplification signal at longer wavelengths. For example, the optimal frequencies 
for detecting the signals shown in Figure Ic and Id are 10^ GHz (sub-mm) and 10^ GHz 
(radio), respectively. The achromaticity of the amplification peak can be used to separate 
the lensing signal from noise due to intrinsic variability or interstellar scintillations. Detailed 
observations of future afterglows are necessary in order to assess the characteristic level of 
intrinsic variability. 

(b) At early times, the temporal profiles of the lensed and unlensed fluxes have the same 
shape but different amplitudes. During this period, the source can still be regarded as 
pointlike and the offset between the lensed and unlensed curves is set by the point source 
magnification factor at a constant b. The unknown value of b could therefore be inferred 
from this asymptotic offset in amplitude between the lensed and unlensed regimes. 

(c) The maximum amplification occurs at the time t^, when the ring crosses the lens, namely 
when Rs ~ b. The otherwise unknown source size Rs can therefore be inferred at the time 
i^. By taking the ratio between the ring size and the period t^,, one finds the mean velocity 
of the expanding ring during that time interval in units of the Einstein radius, te- Given a 
probability distribution for (based on a reasonable mass and redshift distribution for the 
lenses), one could then test the hypothesis of superluminal expansion. 

(d) Analysis of the shape of the lightcurve after the peak can provide more detailed information 
about the fractional width of the ring W and the temporal history of Rs{t). The smaller 
W is, the higher and narrower the amplification peak gets. When Rg ~ b, the value of the 
magnification ^(Rs,W,b) becomes highly sensitive to the source size Rg- By monitoring the 
lensed fiux as a function of time, one could infer the magnification /Uobs(i) = Fohs{t)/ Foit), 
where Fo{t) is found from the power-law extrapolation of the observed -Fobs(i) after the 
end of the microlensing event, to earlier times. Based on the magnification history ^ohs{t) 
one could infer the time evolution of i?s(t) from the constraint fxi^Rg, W, 6) = /^obs 

(t), where 

^i{Rs, W, b) is given by equation (|8|) and b is inferred based on point (b) above. 



The quantitative interpretation of the lensing signatures suffers from an ambiguity about the 
physical size of the Einstein radius of the lens. This ambiguity can be removed through a parallax 
experiment, in which two (or more) telescopes, separated across the solar system, observe the same 
microlensing event with different values of b (Grieger, Kayser, & Refsdal 1986; Gould 1994). Since 



^The variation of the relativistic Doppler effect across the ring might result in a slight chromaticity of the lensing 
signal, but we ignore it here. 
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variability is induced by the superluminal expansion of the source (rather than by the motion of 
the lens, as is usually the case in microlensing events of steady sources), the two telescopes would 
simply observe different lightcurves with different values of b. Based on their known separation 
and their inferred b values [see point (b)] one could then measure the physical size of the Einstein 
radius, te- The shape of the different peaks measured by the two telescopes can be used to test 
for self-similarity in the shock structure and dynamics. 



3.2. Microlensed Polarization 

The different lines in Figure 2 show the deviation from the steady polarization signal that is 



predicted by equation due to microlensing [Eqs. (|Tl| ) and (14)]. The different panels show 
several random realizations for various choices of the lens-source separation b, and the number of 
ring segments A^. We consider two values of N, one in which the ring is composed of a single 
radial strip composed of nearly square segments {N = 63), and a second in which it is divided 
into two such strips (A'^ = 250). The particular value that the polarization obtains at any given 
time t depends on the specific set of random orientation angles {i;^^}^]^ that were assigned to the 
segments in each realization, and so the fluctuations induced by lensing should be analyzed on a 
statistical basis only. 

The main qualitative characteristics of the lensing signal are: 

(a) The polarization changes around t = in coincidence with the flux amplification peak. 
At that time, the polarization fluctuates because as the ring expands, different segments 
approach the lens (and hence the point of maximum amplification) at different times. At 
any given time, the segment which crosses the lens obtains the largest area in the image 
plane and provides the largest contribution to the overall polarization. The fluctuation rate 
increases as the area of each individual segment gets smaller (or as A^ gets larger), because 
smaller segments sweep faster across the lens. 

(b) If the ring is narrower than the Einstein diameter at lens crossing (i.e. Wb ^ 1), then the 
typical fluctuation amplitude, 5 = {{P) / (Pq)) — 1, is roughly independent of A^ (see top 
panels of Fig. 2). In this case, the ring is sliced into a fixed number of "effective" segments, 
each having a length of order the Einstein diameter, so that N^f^ ~ (27r/3s)/(2rE) ~ nb, and 

(c) The fluctuation amplitude decreases with increasing b, because in this limit the highly- 
magnified zone behind the lens amounts to a smaller fraction of the entire ring area. 
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4. Conclusions 

We have shown that microlensing by stars can be used to study the size, superluminal 
expansion rate, and granularity of the photospheres of GRB afterglows. The light curves shown 
in Figure 1 can be used to extract the source impact parameter b relative to the lens (based on 
the normalization offset between the pre- and post-lensing curves), the fractional width of the 
emission ring (from the height and width of the amplification peak), and the source expansion rate 
and size in units of the Einstein radius of the lens. The source size can be measured explicitly 
through a parallax experiment which would obtain two (or more) light curves that sample the 
achromatic amplification peak at different times (cf. Fig. 1). Such an experiment could serve as 
the definitive tool for discriminating between a microlensing event and intrinsic variability of the 
afterglow source. 

By monitoring the variability of the polarization with time during a microlensing event, 
it is also possible to estimate the number of coherent magnetic field patches on the afterglow 
photosphere (Fig. 2). 

If the cosmological density parameter of stellar mass MACHOs is O^, then most afterglow 
events will acquire an impact parameter b ^ 10(r2*/0.1)~^/^ from their nearest lens. Multi-band 
photometry with an accuracy of ~ 0.03 mag, could then detect the flux amplification signal shown 
in Figures la-c and test for its achromaticity, or else place interesting upper limits on O*, based on 
a relatively small sample of frequently-monitored afterglows. The ~ 1% amplification signal shown 
in Figure Id for 6 = 10 would appear in 5-100% of all afterglows after 2-3 months, at the time 
when the peaJi flux of ~ mJy is reached in the radio, at ~ 10^ GHz. A future X-ray satellite which 
would locate afterglows to within an arcminute (like BeppoSAX does) for all GRBs detected by 
BATSE, might identify hundreds of afterglows per year, and could provide a rich sample for such 
microlensing studies. 

Although our results were limited to isolated point lenses, their qualitative features should be 
common to lens systems with more complicated caustic structure, such as binary stars or galactic 
cores. 

We thank Eli Waxman for valuable discussions and for communicating results from his work 
prior to publication, and David Heyrovsky for useful comments on the manuscript. This work was 
supported in part by the NASA ATP grant NAG5-3085 and the Harvard Milton fund. 
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Fig. 1. — The unlensed (solid line) and lensed (broken lines) flux from a GRB afterglow with £'52 = 1 
and ni = 1 at a frequency v = 10^^ Hz. The different broken lines correspond to different fractional 
widths of the emission ring, W = 5% (short-dashed, highest peak), 10% (long-dashed, middle peak) 
and 20% (dot-dashed, lowest peak). The lens mass is M = IMq and its redshift is z/ = 0.5. The 
source redshift is Zs = 2. The likelihood for the events shown is ~ (10-30)%(O*/0.1)(6/3)^ (see 
Fig. 1 in Gould 1995). 
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Fig. 2. — The lensed polarization signal, (P), normalized by the (constant) unlensed value {Pq}- 
The different lines show three random realizations of the time-varying polarization that would be 

observed during a microlensing event if the emission ring is composed of N nearly-square segments 
which produce a polarization of equal amplitude but random orientation. Results arc shown for 
different values of A'' and the source-lens separation b. The unlensed polarization is (Pq) ~ 0.09 
and 0.04 for N = 63 and 250, respectively. Parameters are the same as in Figure 1 with W = 10%. 



